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In this paper we explicate a method of quantum hydrodynamics (QHD) for the study of the quantum 
evolution of a system of polarized particles. Though we focused primarily on the two-dimension physical 
systems, the method is valid for three-dimension and one-dimension systems too. The presented method 
is based upon the Schrodinger equation. Fundamental QHD equations for charged and neutral particles 
were derived from the many-particle microscopic Schrodinger equation. The fact that particles possess 
the electric dipole moment (EDM) was taken into account. The explicated QHD approach was used to 
study dispersion characteristics of various physical systems. We analyzed dispersion of waves in a two- 
dimension (2D) ion and hole gas placed into an external electric field which is orthogonal to the gas plane. 
Elementary excitations in a system of neutral polarized particles were studied for ID, 2D and 3D cases. 
The polarization dynamics in systems of both neutral and charged particles is shown to cause formation of 
a new type of waves as well as changes in the dispersion characteristics of already known waves. We also 
analyzed wave dispersion in 2D exciton systems, in 2D electron-ion plasma and 2D electron-hole plasma. 
Generation of waves in 3D system neutral particles with EDM by means of the beam of electrons and 
neutral polarized particles is investigated. 



I. INTRODUCTION 

A two-dimension electron gas (2DEG) is in the focus of 
many studies. Sometimes it is considered to be an element 
of the spin-field effect transistor [1]. Ions and holes are 
distinct from electrons in that they can polarize in an exter- 
nal electric field or arrange in the line of the external field 
in the case they are rigid dipoles. EDM of ions and holes 
may cause changes in the dispersion of known eigenwaves 
as well as the onset of waves of novel types. 

Dispersion characteristics of 2DEG have been analyzed 
in works of 0-@l- Magnetoplasma waves in 2DEG have 
been studied in 12, (3D- Excitation of charge and spin den- 
sity in 2DEG has been addressed in |0] where spin-orbital 
interaction has been taken into account and the 2DEG has 
been supposed to locate in an external magnetic field. Not 
only wave processes attract researchers' interest, but 2D 
magneto-transport also does. Studies of the effect of spin- 
orbital interaction on the 2D magneto-transport are also im- 
portant. An equation set has been presented in [5] that de- 
scribes charge and spin diffusion in 2DEG with account 
of spin-orbital interaction [6]. 2D hole gas along with 
2DEG is used for construction of transistors, for example 
Atomic Layer Doping-Field Effect Transistor (ALD-FET) 
13] ■ Collective modes in 2D bilayer graphene are consid- 
ered in article [8]. 

In recent years attention is paid to the effect of the in- 
trinsic magnetic moment (IMM) on the characteristics of 
charged particle systems. The propagation of perturbation 
of IMM and EDM does not require much energy as it oc- 
curs without mass transfer. Both processes may be used in 



the information transfer. In biological systems, for exam- 
ple, polarization processes, i.e. EDM propagation, are the 
predominant way of signal transfer. 

There are modern studies which are focused primarily 
on the EDM dynamics 13], [flQII . They deal with quasi- 
2D and multilayered systems of both charged and neutral 
polarized particles. There are articles where take into ac- 
count of the EDM influence on the Bose-Einstein conden- 
sate (BEC) iflXTI , 1121 . In a particular work a contribution 
of polarization to the dispersion of Bogoliubov's mode has 
been analyzed. 

Influence of EDM on dynamic of charged particles is not 
large, analogously with the influence of IMM. But, evo- 
lution of IMM lead to existence of new physical effects 
which described below. In this way we can expect the ex- 
istence new physical process which are the consequence 
of EDM dynamics. IMM dynamics causes changes in the 
wave dispersion in magnetized plasma IU3H15I1 as well as 
the existence of novel branches of dispersion in systems of 
suchkind [14]. Interaction between EVIM of neutron beam 
and magnetized plasma leads to the generation of waves in 
plasma llal. lfl5ll . This effect is due to spin-spin, spin- 
current IU5I1 and spin-orbital interactions [16] between 
IMM of the beam and IMM and currents in plasma. A lot 
of works is also devoted to the dynamics of IMM in BEC. 

Classical methods were used previously to create a de- 
scription of the collective dynamics of charged particles 
that takes into account the EDM that arises in the medium 
as a consequence of charged particles' movement II 1711 . 

A quantum mechanics description for systems of N 
interacting particles is based upon the many-particle 



Schrodinger equation (MPSE) that specifies a wave func- 
tion in a 3N-dimensional configuration space. As wave 
processes, processes of information transfer, diffusion and 
other transport processes occur in the three-dimensional 
physical space, a need arises to turn to a mathematical 
method of physically observable values which are deter- 
mined in a 3D physical space. To do so we should de- 
rive equations those determine dynamics of functions of 
three variables, starting from MPSE. This problem has 
been solved with the creation of a method of many -particle 
quantum hydrodynamics (MPQHD). Therefore, the reason 
for development such method is analogous with the moti- 
vation of the density functional theory [18]. In our work 
we propose a further development of the MPQHD method. 
Here we consider the EDM dynamics in systems of charged 
and neutral particles. The MPQHD method has been ap- 
plied to systems of charged I119l - l2lll . II22I1 an d neutral I123ll 
particles before. Many researchers ltl3ll24ll25ll . ||2611 e.g. 
Madelung [27] and Takabayashi [28] consider derivation 
of QHD equations from a Schrodinger equation for a sin- 
gle particle in an external field. Marklund et. al. ltl3ll 
suggests the method of generalization of one-particle QHD 
equation for description of many-particle systems. In ar- 
ticle IJ26I1 method proposed in ltl3ll is used for derivation 
of QHD equations for system of relativistic particles. Both 
classical and quantum dynamics of a plenty of interacting 
particles in a configuration space are discussed in II29I1 . 

Bibliography contains many examples of QHD applica- 
tion for the analysis of various processes and phenomena. 
Ion-acoustic waves in dusty plasma Il30l432ll and the dis- 
tribution of non-linear electrostatic solitary excitations in 
it 03il have been analyzed. QHD method is also used to 
study wave processes in electron-positron-ion plasma 134 - 
36T1 . A system of QHD equations can be applied to the anal- 
ysis of instability of quantum plasma ll37H4lll . and in par- 
ticular of the modulation instability of electron cyclotron 
waves [38] and magnetoacoustic waves [39]. Dispersion 
characteristics of magnetoacoustic waves have been ana- 
lyzed in this way too I42l - l44ll . The dispersion of elec- 
trostatic waves with frequencies below the electron plasma 
frequency has also been studied using QHD [45]. The dis- 
persion relation for ion acoustic waves and the absorption 
coefficient for Landau damping are obtain in paper D8I1 
and the quantum electrodynamical short wavelength cor- 
rection on plasma wave propagation for a nonrelativistic 
quantum plasma is investigated in article [46]. 

Some of the results that have been obtained using QHD 
equations derived from a single-particle Schrodinger equa- 
tion are presented in the review [47]. 

The MPQHD approach to systems of neutral particles, 
e.g. of ultracold boson-fermion mixtures, has been devel- 
oped in 112311 and applied later to describe linear and non- 
linear characteristics of BEC and of boson-fermion mix- 
tures EiiQ. 

In this work we extend MPQHD approach to a sys- 
tem of EDM-having particles. The method we developed 



here may be an effective tool for the investigation of static 
and dynamic behavior and transport characteristics of in 
graphene [51] and nanofluidics [52], the last one deals 
with the dynamics of ions in water, a substance with large 
EDM of molecules. 



In this article, using the QHD approach we calculate the 
dispersion of polarization waves and spatial charge waves. 
We show the EDM dynamics to cause novel branches of 
dispersion in quasi-2D plasma-like media and in systems 
of neutral particles. We also accomplish analytical calcula- 
tions of dispersion characteristics for the waves we discov- 
ered. 



Electrically polarized particle system can interact with 
the beam of charged and polarized particles by means 
charge-dipole and dipole-dipole interaction. Such interac- 
tion could lead to transfer energy from beam to medium 
and, consequently, to generation of waves. In plasma 
physics the effect of generation of waves by means elec- 
tron 115311 - Ip4II or magnetized neutron [16] beam is well- 
known. In presented article we consider similar effects in 
system neutral polarized particles. 



Our paper is organized as follows. In Sect. 2 we 
present the derivation of the momentum balance equation 
for EDM-having charged particles from MPSE in a self- 
consistent field approximation. An explicit form of the 
quantum part of the pressure tensor is also presented. In 
Sect. 3 we obtain equations of polarization evolution and 
polarization current. The self-consistent field approxima- 
tion is used. In Sect. 4 a calculation is accomplished of 
eigenwaves in a 2D system of EDM-having charged parti- 
cles located in uniform external electric field. New disper- 
sion branch of w(k) is shown to exist and the contribution 
of polarization into the dispersion of 2D Langmuire wave 
is estimated. In Sect. 5 we show the existence of a po- 
larization wave along with acoustic wave in a system of 
neutral polarized particles for ID, 2D and 3D cases. In 
Sect. 6 dispersion characteristics of a two-sort 2D sys- 
tem of charged particles are discussed. An assumption is 
made that particles of one of the sorts bear EDM. We show 
that polarization dynamics here leads to the existence of 
a new dispersion branch. Analytical relation for w(k) is 
constructed. Deep analysis of quantum magnetoacoustic 
waves is performed. In Sect. 7 we apply the results pre- 
sented in sections 4-6 to the analysis of excitations in 2D 
electron plasma and in a system of excitons. In Sect. 8 we 
show that there is instability at interaction of beam of po- 
larized particles with the polarized medium. The increment 
of instabilities is calculated. In Sect. 9 we obtain the incre- 
ment of instabilities which arise at electron beam propaga- 
tion through system of neutral particles with the EDM. In 
sect. 10 brief summary of obtained results are presented. 



II. CONSTRUCTION OF FUNDAMENTAL EQUATIONS 
AND THE MODEL ACCEPTED 

In this section we derive the MPQHD equations from 
MPSE. Here we present the key steps of getting the 
MPQHD equations. We receive the equations for the sys- 
tem of charged particles with EDM. Obtaining equations 
could be used for neutral particles with the EDM as well. 
Method of MPQHG allow to present dynamic of system 
of interacting quantum particles in terms of functions de- 
fined in 3D physical space. It is important at investigation 
of wave process, which take place in 3D physical space. 

Starting from MPSE 



—ihd t tp = Hip 



with the Hamiltonian 
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i D? + ^ 



d a E a 



+ 



i,3=£i 



le^Gv + e^CZ-^G" 



(1) 



we construct a system of QHD equations for charged par- 
ticles that have EDM df. Equations constructed in this 
way are also valid for neutral polarized particles. The 
following designations are used in the Hamiltonian {T): 
Df = -ihdf-e t Af ext /c, <Pi, ext , Af ext - the potentials of 
external electromagnetic field, Ef ext -, is the electric field, 
quantities e,-, m^-are the charge and mass of particles, h-is 
the Planck constant, and Gy = 1/r^, Cf- = —dfl/rij, 

Gff = dfd^l/rij - the Green functions of the Coulomb, 
charge-dipole, dipole-dipole interactions, respectively. 

The first step in the construction of QHD apparatus is 
to determine the concentration of particles in the neighbor- 
hood of r in a physical space. If we define the concentra- 
tion of particles as quantum average of the concentration 
operator in the coordinate representation h = J2i <K r — r i) 
we obtain: 

n(r,t)= f dR^25(r-r^*(R,t)^{R,t) (2) 
■' i 

where dR = Yl p=1 dr p . 

Differentiation of n(r, t) with respect to time and apply- 
ing of the Schrodinger equation with Hamiltonian £Q leads 
to continuity equation 

d t n(r,t)+Vj(r,t) = (3) 

where the current density takes a form of 



j a (r,t) = J d Rj2S(r-r t )^-(^*(R,t)(Dl 



+ {D^{R,t)Y^{R,t) 



4>(R,t)) 



(4) 



The velocity of i-th particle \i(R,t) is determined by 
equation 

\i{R, t) = —V t S{R, t) - —A hext . (5) 

rrii rriiC 

The quantity Vi(i?,t) describe the current of probability 
connected with the motion of i-th particle, in general case 
Vi (R, t) depend on coordinate of all particles of the system 
R, where R is the totality of 3N coordinate of N particles 
of the system R = (r 1; ...,r N ,t). 

The S(R, t) value in the formula (O represents the phase 
of the wave function 

i/)(R,t) = a(R,t)ex P C S ^ R,t " > ). 

h 

Velocity field v(r, t) is the velocity of the local centre of 
mass and determined by equation: 

j(r,t) = n(r,t)v(r,t). (6) 

This means that Ui(r, R, t) = Yi(R,t) — v(r, t) is a 
quantum equivalent of the thermal speed. 

A momentum balance equation can be derived by differ- 
entiating current density d4]l with respect to time: 

8 t j a (r,t) + -^n^(r,t) = F a (r,t), (7) 
m 

where F a (r, t) is a force field and 

n Q/3 (r,t)= f dRJ2 S ( r ~ r i) x 



^ + (R,t)(DfD^)(R,t)+ 



Arrij 



+ (Dfif;) + (R,t)(DfiP)(R,t) + c.c. 



(8) 



represents the momentum current density tensor. 

Let's now perform explicit separation of particles' ther- 
mal movement with velocities U;(r, R, t) and the collective 
movement of particles with velocity v(r, t) in equations of 
continuity (0I and of the momentum balance (JT]). We can 
see now that the tensor Il a/3 (r, t) takes the form 

n a ^(r, t) = mn(r, t)v a (r, t)v^(r, t) 



+p ap {r,t) + T ap {r,t) 



a/3 I 



In this formula 



p afi (r,t)= fdRj2K r - r iy(R,t)miU?u? (9) 

is the tensor of kinetic pressure. This tensor tends to zero 
by letting Uj — >• 0. 
The tensor 



T ap {r,t) 



2m 



r N 

/ dRJ2 s ( r - r i) a2 ( R ,t 



d 2 \na(R,t) 



(10) 



is proportional to H 2 and has a purely quantum origin. For 
the large system of noninteracting particles, this tensor is 



T#(r,t) = -—d a d p n(r,t) 
Am 



n 2 



i 



+ —^^(d a n(r,t))(d p n(r,t)). (11) 

This term is named Bohm quantum potential. 

As the particles of the system under consideration in- 
teract via long-range forces the approximation of a self- 
consistent field is sufficient to analyze collective processes. 
With the use of this approximation two-particle functions 
in the momentum balance equation can be split into a prod- 
uct of single -particle functions. Taken in the approximation 
of self-consistent field, the set of QHD equation, continuity 
equation and momentum balance equation has a form: 

d t n{r,t)+V(n(r,t)y(r,t)) = 0, (12) 

mn(r, t)(d t + vV)u°(r, t) + fy(p a/J (r, t) + T al3 (r, t)) 
= en(r, t)E a ext {r, t) + P*(r, t)d a EU*, t) 
+en(r,t)e a ^v^r,t)B2 xt (r,t) 

-e 2 n(r,t)d a f dr'G(r,r')n(r',t) 
-en(r,t)d a d p f dr'G(r,r')P^(r',t) 
+eP p (r 1 t)d a d f dr'G{r,r')n(r',t) 

+ pP{r,t)d a f dr'G^(r,r')P^{r',t). (13) 

Let's discuss the physical significance of terms on the 
right side of ( fT3l ). The first three terms describe the inter- 
action with external electromagnetic field. The first term 
represents the effect of external electric field on the charge 
density. The second term is the effect of non-uniform ex- 
ternal electric field on the polarization density. The form of 
this term is similar to that of the force field which affects 
the magnetic moment in a magnetic field 



1200, Q28Q. It 

should be noted that the form of this term is distinct from 
the expression that describes the force affecting a single 
dipole. The third term of the equation (PT3l l is the Lorentz 
force. Other terms in (fT3l l describe a force field that repre- 
sents interactions between particles, namely the Coulomb 
interaction of charges, the effect of dipoles on charges, 
charges on dipoles and dipole-dipole interactions. 



Note, that for a 3D system of particles the momentum 
balance equation (fT3l l may be written down in terms of 
electrical intensity of the field that is created by charges 
and dipole moments of the particle system: 

mn(r, t)(d t + vV)/(r, t) + d (p a ^(r, t) + T^(r, t)) 



en(r, t)E a (r, t) + P"(r, t)d a E (r, t) 



+ en(r,t)£ a0 ' r v (r,t)B, 



7 

ex I. 



El 



(14) 



where: E a {r,t) = E« xt (r,t) + E? nt {r,t) and 
£7? t (r,t) = E*(r,t) + E°[(r,t). These variables 
meet equations divE q (r,t) = Anp and divE d (r,t) = 
—4:irdivP(r,t) where p = J2 a e a n a( r ,t). This leads to 
a field equation 



divE int (r, t) = 4np — AtrdivP. 



(15) 



The method we develop in this work is valid both for 
bosons and fermions. The type of statistics that particles 
are subject to affects the calculation of many -particle func- 
tions (correlations) that evolve in the momentum balance 
equation (TT3l and are neglected in the self-consistent field 
approximation. A method for the calculation of correla- 
tions in QHD equations has been developed in works II19I - 

Polarization evolves in the momentum balance equation 
GSarCS): 

P a (r,t)= f dRj28(r-TiW(R,t)d?il>(R,t). (16) 

■' i 

The EDM operator df there affects coordinates of the 
i-th particle. The polarization thus affects the evolution 
of concentration and of the velocity field. To calculate 
P a (j,t) may be interesting itself. It is, therefore neces- 
sary to derive an equation that describes the evolution of 
polarization. 



III. EQUATIONS FOR THE EVOLUTION OF 
POLARIZATION 

To close the QHD equations set (fT2l) . ( [TSl we derive 
equation for the polarization evolution. If we differenti- 
ate the definition for polarization ( TToT) with respect to time 
and apply the Schrodinger equation, the required equation 
for the polarization evolution can be obtained: 

d t P a (r,t)+d R a ?{r,t) = 0. 



(17) 



A polarization current density 

R afs (r,t) = [ dRY^r-r, 



2m, 



x i,*(R, t)(D^(R, t)) + (Dfy(R, t))*i>{R, t] 



(18) 



occurs in this equation. 

We have two ways to close the QHD equations set. The 
first one is to express i? Q,3 (r, t) in terms of n(r, t), v a (r, t) 
and P a (j,t) using additional assumptions or experimen- 
tal data. The other way is to derive the equation for evo- 



lution R a,3 (r, t) in the same fashion it was accomplished 
previously for other material fields, for example quantities 
considered above n(r, t), v a (r, t) and P a (r, t). Now the 
evolution equation R al3 (r, t) occurs in the form of 



8 t R a ^r,t) + -8rrR a ^(r,t) = -P a (r,t)E p ext (r,t) 



m 



m 



e 

m 



+-^s^ s R^(r,t)B 5 ext (r,t) + -(d El ot {r,t)) f di?J>(r - r p )d%(F p P{R,t)il>{R,t) 
mc m J ^-^ p ' 

p 

dr'(d?G(r,r')) J dR ]T S(r - r p )6(r - r n )d%P(R,t)il>(R,t) 
dr'(d G^(r',r)) / dRj2 S(r - r p )5(r - r^cD,^* (R, t)1>(R, t) 

P,Uy^p 

+^-[dr'(d?G->(r,r'))[dR ]T 5(r-T p )6(r-r n )<%W(R,t)il>(R,t) 
-^Jdr'(d^G^(r,r'))JdR ]T <S(r - r p )8(r - rJd^X^CMVK^)- 



e 
m 



(19) 



p.n^p 



Now we consider the physical meaning of the terms in two-particle correlation function. The two-particle func- 
right side of this equation. The first three terms describe tion found in the four last term of equation ( fT9l . In articles 
the interaction of particles with external electromagnetic 
field. The first term represents the effect of external elec- 
tric field on the charge density. The second term is the methods developed in 11190 . Il2lll show that statistics affect 



1911 . Il2lll the method of calculation of correlation is devel- 
oped. Direct calculation of many-particle function using 



effect of non-uniform external electric field on the polar- on correlations only and have no influence on terms arisen 

ization density. The form of this term is similar to that in self-consistent field approximation. Below, in this article 

of the force field which affects the magnetic moment in a we bound our self by self-consistent field approximation, 
magnetic field. The third term of the equation ( fT9l l is an We can see that the equation (T9h contains information 

analog of the Lorentz force. Other terms in (fT9l l describe about the effect of interaction on the transformation of 

a tensor field that represents interactions between particles, R a/3 (r, t) and, as a consequence, on the polarization evo- 

namely the Coulomb interaction of charges, the effect of lution P a (r, t). 



dipoles on charges, charges on dipoles and dipole-dipole 
interactions. 
We used the denotation in (O for R a/3 ^(r, t): 

R a ^{r,t) 



In the approximation of the self-consistent field, the 
equation (fT9l l takes the form 

d t R a ?(r,t) + -dm a ^(r,t) = -P a {r,t)E p ext (r,t) 



m 



m 



dRj2$(.i 



df 

4m.,- 



^{R,t){D p i D^){R,t)+ 



+ (D^) + (R,t)(D^)(R,t)+c.c). (20) 



Equation (fT9l l show the evolution of polarization current 
i? Q/3 (r, t) in consequence of interaction of particles with 
external and internal force field. Equation ( fT9l l contain 



+_L e ^*U«*7( r? t)B s ext (r, t) + -D^(r, t)d? ET ext (r , t) 
mc m 



-—P a (r,tW f dr'G(r,r')n(r',t) 
m J 

D a ^(r,t)dP J dr'G 7 (r',r)n(r',t) 



e 
m 



+— P a (r,t)d? fdr'G J (r,r')P^(r',t) 

7TI J 

+ —D ar '(r,t)d f> f dr'G^(r,r')P 5 (r',t). (21) 
m J 

The physical meaning of terms on the right side of this 
equation is similar to that of terms in the equation ( fT9l ). 
We used the denotation 

D a ?(r,t) = fdRJ2K r - r i)didi^*(R,t)^(R,t) 



in the equation d2"TT) . 
It can be assumed that 



D*(r,t )oi *r<r>W) 



(22) 



(23) 



n(r,t) 

where a- is the nondimension constant. 

Likewise the previous section we can reformulate the 
equation d2~TT l for the three-dimension case in terms of the 
electrical field ( fT31) . The equation (|2~TT | in this situation has 
the form 

d t R a0 (r,t) + —d^R a ^(r,t) = —P a (r,t)E ? (r,t) 
m m 

+ —£^ s R a ^(r, t)B 5 Ar, t) + — D a7 (r, t)d p E^(r, t). 
mc m 

(24) 

If we put a velocity field into the equation d2~TT l the tensor 

i? Q ° 7 (r, t) transforms into 

R a ^(r, t) = r Q/37 (r, t) + T a ^(r, t) 
+mR afS (r, t)u 7 (r, t) + miT 7 (r, t)v p (r, t) 



mP a (r,t)v (r,t)v~<(r,t). 



(25) 



In this representation the contribution of thermal move- 
ment 



„«/37 



N 

r,t)= / dR^2S(r-r t )d^a a (R,t) 



8 7 



(26) 

becomes explicit and the analog for the Bohm quantum po- 
tential occurs 



fi 2 r N 



r - r, x 



xdy\R,t) 



d 2 In < 



(27) 



dxpidx^ ' 

Approximate connection of T Q/37 (r, t) with the concentra- 
tion and polarization of particles is presented in work II 1211 . 
The major contribution into alterations of polarization in 
a system of charged particles is from charge interactions 



and from the effect that the external electrical field has on 
charges. As a result 



d t R^(r,t) = -P a (r,t)E p ext (r,t) 
m 

+ —P a (r, tW [ dr'G(r, r>(r', t) (28) 
m J 

can be derived from the equation d2~TI) . 

Later in this paper the equation ( |28l is used to analyze 
the elementary excitations spectrum in a two-dimensional 
system of EDM-having charged particles. 

Method of QHD allows derive an one-particle nonlinear 
Schrodinger equation (NLSE). Corresponding one-particle 
wave function describes evolution of many-particle sys- 
tem in 3D physical space. Most famous NLSE are the 
Ginzburg-Landau equation in superconductivity 15511 and 
the Gross-Pitaevskii equation in the physics of ultracold 
gases and Bose-Einstein condensate [56]. In this article 
we consider the derivation of NLSE. We do not use NLSE 
for studying of wave dispersion and wave generation. Be- 
cause of it we present NLSE in Appendix. In this paper we 
receive NLSE for the charged particle having EDM. Inter- 
esting consequence of derivation of NLSE is the obtaining 
of London equation. Derivation of the London's equation 
from QHD also is presented in appendix 



IV. THE EIGENWAVES IN A 2D SYSTEM OF CHARGED 
PARTICLES WITH EDM 

One of the primary goals of this article is to derive dis- 
persion characteristics of eigenwaves in 2D systems of neu- 
tral and charged particles that take account of the EDM col- 
lective dynamics. In this section we consider 2D systems 
of charged particles. To do that let's analyze small pertur- 
bations of physical variables from the stationary state. The 
ionic or holes gas is assumed here to localize on a x, y- 
plane 

n = n + Sn, v a = + 5v a , P a = P Q + 5P a , 

(29) 

K« = E o$ az , Po = kE S = ^o£ az , 
where n is a surface concentration. 

We can mark, in this section and below, all physical 
quantity is presented in the form of sum of equilibrium part 
and small perturbations 

f = fo + Sf. 

In this case if we assume that linear excitations 5f are 
proportional to exp{—iujt + zkr) a linearized set of equa- 
tions (fT2l . ( fTSl . (flTt . d28l) gives us the dispersion equation 
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FIG. 1. The figure presents the dependence of the variable f3(k) 
on the wave vector k. Ionic or atomic radius tq is assumed to 
equal 0. 1 nm. 



where 



P(k) = 2vr 



dr 



Mr) 



here £ = r k,r there is an ionic or atomic radius and k = 

^T^ is a modulus of the wave vector. As X mm = 

2ir/k max > 2r then £ C (0, %). The function /3(£) is, 
therefore, positive. Its explicit form is presented at Fig. 1. 
Solutions of the equation d30l) are presented by two 
branches of the dispersion characteristic 
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FIG. 2. The figure shows the dispersion characteristic of the 
Langmuir wave frequency u versus the wave vector k, which 
is described by the equation 03l >. The ionic radius vq is sup- 
posed to be 0.1 nm. Equilibrium polarization has form Pq = 
kEq. Static electric permeability k is defined by the equation 
k = nop 2 /(3kgT). po - is a dipole moment of an ion or atom, 
T - temperature of the medium, kg - Boltzmann constant. Sys- 
tem parameters are assumed to be as follows: no = W s sm~ 
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(31) Po = 3-l(T^Cm,T = 100K, E = 3- KTV/m, m l = W~^g 
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(32) 
where w Li2 = 2Tre 2 n k/m denotes the Langmuir 2D fre- 
quency. 

Provided that polarization effects are small the square 
root expression in (l32l may be expanded to yield formulae 
co(k) in the following form 
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The relation (1331) expresses the dispersion of 2D Lang- 
muir waves with the account of dipole-dipole interactions. 
The solution (l34l expresses the dispersion of waves that 
emerge as a result of polarization dynamics. Plots at Figs. 
2 and 3 represent relationships (l33l . d34l) . respectively. 



FIG. 3. The figure shows the dispersion characteristic of the 
quantum ionic polarization wave frequency cj versus the wave 
vector k, which is described by the equation < !34b . The ionic ra- 
dius ro is supposed to be 0. 1 nm. Physical conditions of the sys- 
tem are supposed to be the same as described in the legend to Fig. 
2. 



V. EIGENWAVES IN A SYSTEM OF POLARIZED 
NEUTRAL PARTICLES 

In this section we consider ID, 2D and 3D systems of 
EDM-having neutral particles. To do that we use the equa- 
tion ( f2Tb . Since the system of neutral particles resides in 
a uniform electromagnetic field, there is only one term on 
the right side of the equation (|2~TT) . It is also assumed that 
interactions make the largest contribution into the changes 
inR afi (r,t). 

If so then the equation d2ll and (l23l transforms into 

d t K p {r,t) = a — x 

mn(r, t) 
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[ dr'G^i 



r,r')P s (r',t). 



(35) 
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FIG. 4. The figure shows the dispersion characteristic of the 
2D quantum atomic polarization wave frequency u>(k) versus the 
wave vector k, which is described by the equation d36l l. The 
atomic radius is supposed to be 0.1 nm. System parameters are 
the same as in Fig. 2. 



If taken in a linear approximation, the set of equations 
([T2l l. (O, ( fTTT l. ( 1351 ) for neutral particles is closed. This 
allows the analysis of polarization waves in a system of 
neutral particles. Equations of continuity (fl2l l and of the 
momentum balance (ITS! herein describe the dynamics of 
acoustic wave. If we derive a solution for eigen waves in a 
2D system the dispersion equation has a form of 
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where /3(k) is defined by the relation (l3~T1 l. The dispersion 
dependence (l36l l is presented on Fig. 4. 
In ID case u>(k) occurs as 
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where 
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The quantity (l38l l and ID dispersion dependence (l37l l are 
presented on Figs. (5) and (6) correspondingly. 

The set of equations ( fTTl l, d35l , ( f23l may be applied to 
analyze 3D systems, otherwise the electric field E(r, i) that 
dipoles create may be introduced explicitly ( fT5t , ( fTTl l, (f24l . 
(fH. 

As this is done the dispersion equation to(k) transforms 
into 
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Equations (l36l l and (l37l) differ in the power coefficient 
by the wave vector k and also with the coefficients f3{k) 
and /?! (/c), which depend on the wave vector, occur in ID 
and 2D cases, respectively. 



FIG. 5. The figure presents of the variable /3i(fc) on the wave 
vector k. Ionic or atomic radius tq are assumed to equal 0.1 nm. 
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FIG. 6. The dependence of frequency u> on the wave vector k 
is displayed for the case of single dimension polarization mode 
which dispersion characteristic is defined by the equation < !3Tb . 
The atomic radius ro is supposed to be 0.1 nm. Equilibrium 
polarization has form Po = k-Eq. Static electric permeabil- 
ity k is defined by the equation n = noPQ/(3fcsT). po - is 
a dipole moment of an atom, T - temperature of the medium, 
ks - Boltzmann constant. System parameters are assumed to be 
as follows: n = 10 4 sm~ 1 , p = 3 ■ 10~ 20 Cm, T = WOK, 
E Q = 3- 10 4 V/m, m = W~ 23 g. 



VI. EIGENWAVES IN A 2D TWO-SORT SYSTEM OF 
CHARGED PARTICLES 

A 2D system comprised by charged particles of two sorts 
can be modeled either in a thin film of metal and semicon- 
ductors or in a thin, nano-sized ionic crystal. Only ions and 
holes can have EDM in the former case. Both sorts of ions 
can have EDM in the case of ionic crystal. In our analysis 
we deal with a material where one sort of ions has much 
greater EDM than the other and investigate elementary ex- 
citations in such system. 

The QHD equation set for the system in question is com- 
prised by continuity equations (fT2l) and momentum balance 
equations © for both sorts of particles and also by equa- 
tions (fTTJ, (fJJJ , which describe the polarization dynamics 
of particles those having greater EDM. The field on the 



right side of equations © and d2TT) comes from particles of 
both sorts. Likewise previous sections the 2D system is as- 
sumed here to be located in an external electric field which 
is orthogonal to the plane of particle movement. Solution 
of the eigenwave problem in a two-sort system of charged 
particles leads to the dispersion equation 
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(40) 



If we analyze a film of metal the indices e and i relate to 
electrons and ions, respectively. For the case of semicon- 
ductors index i describes a holes. In the case of a nanolayer 
ionic crystal i relates to ions of different sorts. It is assumed 
in both cases that particles denoted by the indices e and i 
are those having EDM. 

If there is no polarization then the equation (1401 is a 
second-degree equation with respect to lo 2 . Its solutions 
have the form 
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(41) 



The following denotation is used here: 

v 2 qsa =v 2 a + h 2 k 2 /Am 2 a , 

where a denotes e for electrons and i for ions. 

The equation d4TI) contains dispersion characteristics for 
the Langmuir wave and ion-acoustic wave. In the quantum 
case we consider these characteristics take on the form of 
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formula (l46l take form: 
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In this approximation exist simple relation for ion-sound 
waves (l43l. 



The solution (1431 differs from the classical case in the 
presence of a term proportional to H 2 which occurs due to 
Bohm quantum potential (fTTTl . Figure 7 shows dispersion 
characteristics of ion-acoustic waves in classical and quan- 
tum cases. 

If polarization is taken into account in the equation (1401 
then new solution occurs. Assuming the contribution of 
polarization to be small we may obtain the following dis- 
persion characteristic for this novel solution: 
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If the approximation d44l) is used then the equation 



takes the form 
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FIG. 7. This figure displays the function u>(k) for two- 
dimensional ion-acoustic waves. The upper curve shows the 
dispersion characteristic of quantum ion-acoustic waves which 
are described by equation (PJ5} while the lower curve shows the 
dispersion characteristic of classic ion-acoustic waves which oc- 
cur from 045 l l if the term proportional to h 2 is neglected. Sys- 
tem parameters are assumed to be as follows: no 
m = 1(T 2 V 
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The equation (l45l l for ion-acoustic waves occurs in the 
same approximation. 

The polarization also causes additional contributions in 
solutions (|4J1 to occur. Thus, a generalization of solutions 
(l4J~b that takes account of polarization appears in the fol- 
lowing form: 
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VII. THE EFFECT OF POLARIZATION ON THE 

DISPERSION CHARACTERISTICS OF 2D 

SEMICONDUCTORS 



that acoustic-like waves and polarization waves may be ex- 
cited in exciton systems and their dispersion is described 
by formulaes d36]l, |37]) and (l39l 

In the case electrons and holes move in a semiconduc- 
tor without significant formation of bound state excitations 
in a semiconductor may be regarded to as excitations in a 
two-sort system of charged particles. The term "hole" re- 
gards to ions in the semiconductor, which may have EDM. 
This means that holes may be dealt with as EDM-having 
particles. 

Using the model described above we can apply formu- 
lae from the Section VII to describe the dispersion of ele- 
mentary excitations in 2D electron-hole plasma. Index i in 
formulas (l40l) -(l48l) regards to holes. 



VHI. EXCITATION OF POLARIZATION WAVES BY 
BEAM OF NEUTRAL POLARIZED PARTICLES 

We do not meet experimental works there were obtained 
the polarization waves received in our article. On exam- 
ple of 3D sample we show the way to obtain such waves 
in experiment. A possible way of generation of waves of 
polarization it is the propagation of beam of electrons or 
neutral polarized particles through the sample. Instead of 
the beam can be considered a current of particles through 
the sample, in context of next section it can be the elec- 
trical current. In this and next sections we show there are 
instabilities which led to generation of waves in the sam- 
ple. The effect of waves generation by beam of electrons is 
well in plasma physics 115 311 . 115411 . 

In this section we use equations (Tf2l . ( fT4l . ( TP71) and (l24t 
for each sorts of particles and equation of field ( Tf5l . In 
right side of equation (TOT l the polarization P is a sum of 
polarization of medium P d and polarization of the beam P fc 

P = Pd + Pfc 

The equilibrium state of system is characterized by fol- 
lowing values of the medium parameters: 

n d = n Qd + 5n d , v^ = + v°[, 



The simplest model of semiconductors is a system com- 
prised by particles of three sorts, namely electrons, ions 
and neutral atoms. Atoms and ions in the semiconductor 
may bear significant EDM. Dynamics of semiconductors 
is usually dealt with as dynamics of electrons and holes. 
Electrons and holes can also form bound states - the exci- 
tons. As excitons have EDM, the dipole-dipole interaction 
plays a substantial role in a system of excitons. 

The approach we present in this paper allows us to de- 
scribe the interaction of excitons via their EDMs and inter- 
action of their EDM with an external electric field. Since 
excitons are neutral EDM-having quasi-particles, the result 
we obtained in Section VI may be used to describe elemen- 
tary excitations in such system. We can therefore suggest 






P o a d + 5P«,Rf = + 5Rf 
and values of the beam parameters: 

n b = n ob + 5n b , v% = U5 za + Svg, 
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The polarization P^ is proportional to external elec- 
tric field En- We consider the case then E n = 



[E o sirup,0, E cosf]. In this case the tensor R^ has 
only two unequal to zero elements: Rq% = R^sirup and 
RqI = R ob cos(p. And we consider the small perturbations 
of described here equilibrium state. 
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In the absence of medium we obtain a dispersion depen- 
dence of polarized beam modes 
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At the presence of medium, far from resonant condition 

k z U ~ u>Db, we obtain 
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For the case uj d ~^> k z U from ( |53l we derive 
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Further we consider resonant interaction neutral polarized 
beam with medium. In this situation we present frequency 
in the form 
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k z U + r). 



At resonant condition uj d ~k z U and rj 3> — we have 
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where £ presented by formula d66l l. In the limit 77 <C -^- 
we obtain a formula which is analogous to (l67l l. But in this 
case we need to do some changes, i.e. m e —> m b and oj Le 
to UJ Db . 



IX. EXCITATION OF POLARIZATION WAVES BY 
ELECTRON BEAM 

In previous section we consider the effect of generation 
of the polarization waves in 3D system of neutral particles 
with the EDM. Here we analyse an analogous effect. It 
is a generation of waves by means of a beam of neutral 
polarized particles. 

Let's analyze the interaction of a single velocity beam 
of electrons that moves along z-axis with infinite three- 
dimension medium comprised of neutral particles. Parti- 
cles of the medium are in equilibrium state and polarized 
without any external electric field. This may occur if the 
electron beam moves through a crystal of piezoelectric or 
through a sample of ferroelectric that has residual polar- 
ization. The result of such analysis would also represent a 
physical mechanism of the interaction of an electron beam 
with spatially confined or low-dimension polarized sys- 
tems. The most significant types of interaction in such sys- 
tems are the Coulomb interaction of electrons in the beam, 



the dipole-dipole interaction in the medium and the charge - 
dipole interaction of dipoles in the medium with the elec- 
tron beam. Let's assume that the beam moves alongside the 
polarization vector and the beam's velocity doesn't change 
significantly while it moves through the sample. 

In this section we use equations (fT2l . (fT4l . ( fT71) and (l24l 
for each sorts of particles and equation of field ([131) . 

If so, solution of the QHD equation set leads to the dis- 
persion equation 
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This equation has three solutions. If an account is 
taken of Bohm quantum potential and no electron beam is 
present, solution 

(56) 
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occurs for the polarization wave, while two beam modes 
occur in the absence of the medium 



u = k z U±Ju; Le + — 



(57) 



If the contribution from the Bohm quantum potential is 
neglected and the limit of a low-density beam is assumed 
the solution for beam mode looks like 



U) OikJJ- 
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If we consider the effect of the dipole-comprised 
medium on the beam mode the following solution is ob- 
tained 



UJ ~ kM± 
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This solution is complex and leads to the instability if 
Dj) > k z U i.e. in the case of long wave approximation. 
If UJ d 3> k z U then solution transforms to 



uj = kJJ\ 1 ± 1 



ULe 
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(60) 



Solution ( |59l is valid for a beam moving in a 3D system 
of dipoles if 

k z U / uj d , (61) 

i.e. if there is no resonance between beam mode and polar- 
ization mode. 

Let's see how a beam affects the dispersion of polariza- 
tion wave in the absence of resonance with beam mode. 
The dependence of frequency on the wave vector has the 
form 
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To analyze a resonance interaction that occurs between 
the electron beam and a system of electric dipoles, i.e. the 
situation of ( 162b we should seek for a solution that looks 
like 



uj = k z ll + 7] 
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4m 2 , 
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In this instance we obtain the following equation of a 
frequency shift m 

hk 2 , 



r] ± rf 



hk 2 



2 9 



2Jco 2 + h2k4 



0. 



(64) 



J D ^ Am 2 , 



If T) 3> — then it has the following solution: 
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In the opposite limit of ?] <C 
occurs as 
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So we have shown above in this section that the elec- 
tron beam that moves in a system of dipoles polarized by 
external electric field induces instability and increases am- 
plitude of the polarization mode. It is assumed that the 
electron beam moves alongside of the external electric field 
and the medium is infinite and homogeneous. It should be 
noted that while the electron beam accelerates infinitely in 
such conditions we assumed changes of the beam veloc- 
ity to be small. This was done because our goal was to 
approximate a system of finite length. 



X. CONCLUSIONS 

In this paper we analyzed wave excitations caused by the 
EDM dynamics in systems of charged and neutral particles. 
Method of QHD was developed for EDM-having particles. 
QHD equations are a consequence of MPSE in which parti- 
cles' interaction is directly taken into account. In our work 
we consider the Coulomb, charge-dipole and dipole-dipole 
interactions. The system of QHD equations we constructed 
is comprised by equations of continuity, of the momentum 
balance, of the polarization evolution and of the polariza- 
tion current. In our studies of wave processed we used a 
self-consistent field approximation of the QHD equations. 

Using QHD equations we analyzed elementary excita- 
tions in various physical systems in a linear approximation. 



Waves in a 2D gas of EDM-having charged particles and 
in a gas of neutral EDM-having particles in various physi- 
cal dimensions were considered. A two-sort 2D system of 
charged particles where particles of one sort are assumed 
to have EDM was also analyzed. 

Dispersion branches of a novel type that occurs due to 
polarization dynamics were discovered in those physical 
systems. Furthermore, the effect of polarization on the 
dispersion characteristics of already known wave modes 
was studied. The contribution of Bohm quantum potential 
which is of purely quantum origin was taken into account 
in the calculations. Waves of the electric polarization we 
discovered possess the following feature: their frequencies 
U) tends to zero provided that k — > 0. 

The polarization modes that have been derived in our 
work should be taken into account in calculations of the 
thermal capacity and non-linear dynamics of acoustic and 
ion-acoustic waves. Polarization modes may also con- 
tribute together with phonons to the process of formation 
of Cooper pairs in superconductors. These polarization 
modes may be applied in the construction of devices that 
implement information transfer processes. Excitements in 
polarization may be used as an alternative to spin waves or 
together with them. It also should be taken into account 
that transfer of polarization disturbances plays a major role 
in the information transfer in biological systems. Such 
processes do not require particles of the medium to pos- 
sess EDM as the dynamics of a system of charged particles 
leads to collective polarization P a (r, t). 

The effect of polarization dynamics on the dispersion 
characteristics of exciton systems and of electron-hole 
plasma is discussed in section VIII. The polarization mode 
occurred in those systems too and its dispersion has been 
analyzed. The existence of the polarization mode may af- 
fect characteristics of spintronic and nano-electronic de- 
vices. 

We show the possibility of the process of waves gener- 
ation in the system of polarized particles by means of the 
monoenergetic beam of neutral polarized particles and the 
monoenergetic beam of electrons. 

Hence, in this work we present the advancement of the 
QHD approach to systems of polarized particles and use 
it to show the occurrence of polarization waves in various 
physical systems comprised by charged and neutral parti- 
cles. 



APPENDIX 

Derivation of the non-linear Schrodinger equation 

In this section we discuss the derivation of NLSE for 
systems of charged particles. The derivation is generally 
similar to the derivation of the GP equation for systems of 
neutral bosons in the state of BEC [23]. The NLSE comes 
about from the continuity equation and the Cauchy integral 
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of the momentum balance equation. The Cauchy integral 
exists provided that the velocity field has the form 

v a (r, t) = —d a 6(r, t) - — A a (r, t) (68) 
m mc 

where 6(r, t) is a velocity field potential. 

Starting from QHD equations we can derive an equation 
for evolution of the model function defined in terms of hy- 
drodynamic variables. Thus a macroscopic single-particle 
wave function may be defined as 



$(r, t) = ^n(r, t) exp(-m#(r, t)). (69) 

If we differentiate this function with respect to time and 
apply QHD equations then given the absence of EDM the 
following equation is obtained: 

ihdMr,t)= |-A-|v-^A eit (r,i)) +ep(r,t) 



c' I dr'^Ar- +r(r,t)U(r,t), (70) 



where 
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The equation obtained d70l has the form of NLSE. 

NLSE d70l) describes collective characteristics in a sys- 
tem of many charged particles. This follows from the 
derivation of NLSE and from the definition of a many- 
particle wave function $(r, t) d69l . A nonlinear term 
that is proportional to e 2 depicts the Coulomb interaction. 
Function T(r, t) (fTTT l is the contribution of the kinetic pres- 
sure and does not contain any interaction. 

Now let's write down the NLSE that arise from QHD 
equations (TT2] | and (fT3l) taking into account the EDM of 
particles: 
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here we use designation T(r, t) described by formula (ITO) . 
Also in (l72l exist action of electric field on polarization 
of particles. This effect is presented by function C(r, t), 
which determined by following formula 
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dr'— 

r ^( r '^) 



+ P p (r',t)V fdr"G^(r',r")P^(r",t)\. (73) 

If we deal with a 3D system of particles we can introduce 
a self-consistent electric field into NLSE as that was done 
in equations ( TT2l) . dT4l) and ( fT3t - In such case equations 
d72l ) and ( [731 transform into 

ihdMr,t)= l-^-[V-^A ext (r,t)) (74) 



+&p(T,t)+T(r,t) + C(T,t)\$(r,t). 

In this case C(r, t) has form 
C(r,t) = f r dr'^—P^r',t)VE^r',t), (75) 

E = — V(/? in equations (l74l . d75l) and this field satisfies 
the equation ( fT31) . 

Given the additivity of electromagnetic field the equa- 
tions (1741 . d75l) can be generalized and A ext substituted 
with a vector potential A that is created by both external 
and internal sources. QHD equations ( fl"3l) . (fT4"l) . (L2TT) and 
d24l) can be generalized in the same way. 



Some words about London equation 

In many books the London equation is derived from 
equation which describe the dynamics of one particle in 
external field. But, the London equation describe the prop- 
erties of many-particle system, i.e. system of electrons in 
superconductors. In this section we present a new way of 
derivation of the London equation from QHD equations, 
which describe collective dynamics of particles. Also, we 
obtain generalization of the London equation 

Given the absence of a magnetic field the condition (l68l 
becomes a condition of vortex-free flow. Assuming the ap- 
proximation of a homogeneous medium n = n = const 
an equation 

CTl 

rot(n Q \(r,t)) = °-rot(A(r,t)) (76) 

mc 

can be derived from |68l . 

If we apply equations n v = j, rotA = B, rotB = 

ine/cj, divB = a well-known London equation 



AB(r,t) 
can be obtained where 



4ire 2 n 
— B(r,t) 

mr 



(77) 
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cji e is the Langmuire frequency and A^ is the Londons' 
constant 11571 15811 - 

In the case of an inhomogeneous medium the following 
relationship appears from the equation I 



rotj(j,t) = rot(n(r,t)A(r,t)) 

mc 



■— n(r,t)B(r,t) - — [Vn(r, t), A(r, £)]. (78) 
mc mc 



This gives us a generalization for the equation (177 



AB(r,i) 



4.7re 2 n(r, t) 



mc 



B(r,t) + — [Vn(r,t),A(r,t)]. 
mc 

(79) 
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